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Abstract. The incorporation of two- and three-dimensional 5- function pertur- 
bations into the path-integral formalism is discussed. In contrast to the one- 
dimensional case, a regularization procedure is needed due to the divergence of 
the Green-function G'^^''(x, y; E), (x, y E IR^,IR'^) for x = y, corresponding to 
a potential problem Vipi). The known procedure to define proper self-adjoint 
extensions for Hamiltonians with deficiency indices can be used to regularize 
the path integral, giving a perturbative approach for 5-function perturbations 
in two and three dimensions in the context of path integrals. Several examples 
illustrate the formalism. 



Address from August 1993: Il.Institut fiir Theoretische Physik, Universitat Hamburg, Luruper Chaussee 
149, 22761 Hamburg, Germany. 



1. Introduction 

5-function perturbations play an important role in quantum physics because they can 
serve as simple models for more complicated interactions and, more important, they 
are in many cases explicitly solvable [3]. They are also called point interactions (e.g. 
the Fermi point interaction [9]) in nuclear physics, solid state physics, and for the 
other fundamental interactions. The incorporation of a 5-function interaction in one- 
dimensional models causes no serious problems, and is usually discussed as a simple 
solvable quantum model in many textbooks in quantum mechanics (e.g. [53]). 

The construction of the propagator can be done by a perturbative approach in the 
path integral. Let W{x) = V{x) —^yS^x — a) (a, x G IR), and we consider the 5-function 
as a perturbation of the potential V. Then 

G^^^ {x", x'; E) = G'(^) {x", x'; E) + {T^^} {E)) ~'g^^^ {x", a; E)G^^^ (a, x';E) (1.1) 

with T^l = 1/7 - G^^\a,a]E), and G^^\E) and G^^\E) are the Green functions 
(resolvent kernel) of the unperturbed, respectively perturbed problem. The general 
structure of (1.1) is also known as Krein's formula [3]. 

A comprehensive approach can be found in [32], including numerous examples, 
based on earlier work [7, 27, 51], and, of course, on the monograph of Albeverio et 
al. [3]. Considering the limit 7 —00, i.e. making the strength of the 5- function 
perturbation infinitely repulsive, gives a boundary value problem with Dirichlet 
boundary- conditions at the boundary x = a, therefore explicitly incorporating 
boundary problems in the path integral [15, 16, 18, 36, 37]. Note that for a symmetrical 
model, i.e. V{x) = V{—x), a 5-function perturbation at a; = yields in the limit 
7 — > 00 a doubly generated energy level spectrum with one "spurious ground state" 
\l/o with infinite negative energy and a corresponding wave-function concentrated at 
X = 0, i.e. |\E'o(3^)P = '^(2;) (e-g- [4] and references therein). 

One-dimensional examples, however, describe only quantum mechanical models on 
a line, respectively models which can be reduced to one-dimensional models, e.g. radial 
problems. More realistic ones must be at least two-, respectively three-dimensional. 
They can represent lattices of interaction centers in thin films (e.g. [57]) and solids, 
respectively, and can be chosen arbitrarily located with arbitrary strength. 

Historically, point interactions were introduced in nuclear physics starting in the 
1930s. But this type of interactions appear also in many other branches of physics, 
electromagnetic interactions with vector potential, etc. Also very important, they 
model lattices in solid state physics, yielding in the limit of infinitely many interaction 
centers an electronic band structure. Often the notion of a 5-function perturbation 
is replaced by the notion of particular boundary-conditions which the 5-function 
perturbations actually describe, c.f. [3] for a comprehensive bibliography on this 
subject. 

The purpose of this paper is to incorporate such models into the path integral 
formalism. As we will see, a regularization procedure will be needed. The origin of 
the emerging divergence can be easily seen. A perturbation expansion consideration 
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yields an expression, where the Green function G^^^ (E) of the unperturbed problem, 
say for the free particle Green function, must be evaluated at both arguments being 
equal, i.e. "G''-^''(x, x; i?)". In two dimensions Green functions generally diverge 
logarithmically in this case, i.e. limx^y G(x, y; i?) oc ln|x — y|, whereas in three 
dimensions there is a simple pole, i.e. limx^y G(x, y; oc |x — yl"-*^. Consequently 
a regularization procedure is necessary. It consists basically on a proper definition of 
the corresponding Friedrich extension, making the perturbed Hamiltonian self-adjoint. 
What has to be done is the following: "In physical terms, the coupling constant A in 
the heuristic expression —A + X6y has to be 'renormalized' and turns out to be of 
the form X = rj + arj^, with rj infinitesimal and a e (—00,00]. [3, p. 3]", a topic 
which was first discussed by Berezin and Faddeev [8]. To put it into the language 
of approximating, the (5-function perturbation is rewritten in terms of a Gaussian 
packet ip with parameter e, i.e. ipe{x) = e"*^^"**^ /^/-^Tre, and for D = 3, say, A in 
the expression X'i/j^{x)'ilj^{y)'ilj^{z) has to be explicitly dependent on e, i.e. A = A(e). 
However, this is only necessary for D = 2,3. For D = 1 no regularization is needed, 
and for D > 4, the corresponding Hamiltonian is already self-adjoint [58] , c.f. the note 
in the summary. 

This kind of regularization procedure is known from functional analysis [3] to define 
the corresponding proper self-adjoint Hamiltonian, and our purpose is to show that it 
can also be applied to regularize the corresponding path integral formulation. We will 
find that the entire Green function has to be taken into account in the regularization 
procedure, and not only part of it (compare [50]). 

Such regularization procedures are well-known in other branches of mathematical 
physics, let alone high energy physics, respectively renormalization and quantum field 
theory. For instance, in quantum mechanics let us note the Selberg trace formula [35, 
41, 62, 64], where for non-compact Riemann surfaces, the trace of a properly chosen 
operator-valued function for the Laplacian on the non-compact Riemann surface 
must be regulated by a subtraction of Eisenstein series representing the continuous 
spectrum. This regularization does not simply cancels the continuous part of the trace 
and naively gives zero-terms in the trace formula, but instead leads to inadmissible 
contributions, important in the theory of the Selberg trace formula, and related fields, 
e.g. in number theory. 

The further content of this paper will be as follows. In the next section the 
theory of incorporating two- and three-dimensional 5-function perturbations into the 
path integral will be presented. The regularization procedure will be applied for the 
perturbation expansion. In the third section, several examples illustrate the formalism, 
and will include the free particle, the harmonic oscillator, and the Coulomb potential 
(1/r-) case on the one hand. On the other, some models for nuclear physics potentials 
will be given. The fourth sections concludes with a summary and a discussion of the 
results. In appendix 1 the computation of the Green function for the one-dimensional 
harmonic oscillator is sketched, and in appendix 2 the evaluation of the propagator 
for a free particle with a point interaction in two dimension will be given. 
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2. Time ordered perturbation expansion of the path integral and pertur- 
bation summation for 5-function potentials 



The general method for the time-ordered perturbation expansion is simple. We assume 
that we have a potential W{x.) — V{x.) + V{x.) (x G IR ) in the path integral, where 
it is assumed that W is so complicated that a direct path integration is not possible. 
However, the path integral corresponding to l^(x) is assumed to be known, which 
we call K^^\T). We expand the path integral containing V{x) in a perturbation 
expansion about V^(x) in the following way. The initial kernel corresponding to V 
propagates in e-time unperturbed, then it is interacting with V, propagates again in 
another e-time unperturbed, a.s.o, up to the final state. Let us denote the path integral 
(Feynman kernel, respectively) corresponding to the potential V hj {T = t" — t') 



x(t")=x" 



K(^)(x",x';T)= j P^x(t)expUj^^ 



F(x) 



dt 



x(t')=x' 

We introduce the (energy-dependent) Green function (resolvent kernel) 



(2.1) 



G(^)(x",x';£;) = ^ / dTe^^^/^K(^)(x",x';T) 



K(^)(x",x';T) = -L / G(^)(x",x';i?)e- 



lETh 



dE 



(2.2) 



This gives the series expansion (see also e.g. [7, 23, 27, 32, 51, 61]) 



K(x",x';T) 

x(t")=x" 



x(t')=x' ^ 

K(^)(x",x';T) + 5^( - 



- F(x) - y(x) 



n=l 



1 \ 1 



,n/ ^""-^ 



dt 



dt-i 



X K(^)(xi,x';ti)F(xi)K(^)(x2,Xi;t2 - h) x 



xF(a;„_i)K(^)( ■,tn-tn-l)V{-Kn)K^^\-x", 



n=l 



= K(^)(x",x';T) + E ( - 1 ) n dt, ij-^, 



n n 



X 



K(^) (xi, x'; ti - t')V{xi)K^''^ (X2, xi; t2 - ti) 



X 
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■■■X y(xn_i)K(^)(x,,x„_i;t„-t,_i)F(x,)K(^)(x",Xn;t" -t,) . 

(2.3) 

In the second step we have ordered the time as t' = to < < ^2 < • • • < tn+i — t" and 
paid attention to the fact that K^'^\tj —tj-i) is different from zero only if > tj-i- 
We consider now an arbitrary potential V{x.) in one dimension with an additional 
5-function perturbation so that [32] 

W{x) = V(x) - ^S(x - a) . (2.4) 

The path integral for this potential problem reads 

x{t")=x" 



K^^\x",x';T)= j Vx{t)exp\j^ 



x{t')=x' 



dt) . (2.5) 



Introducing the Green function G{E) of the perturbed system similarly to (2.2), and, 
for the time being, we assume that G^^) (a, a; E) exists. We then obtain due to the 
convolution theorem of the Laplace-Fourier transformation 

G^^\x",x';E) = G^''\x",x';E) ^ ' ' ^ K , , J 

Giy){a,a;E)-l/-f 

Let us consider 5-function perturbations in two and three dimensions. A formal 
series summation produces a formula similar to (2.6). However, a simple example, 
say the free particle Green function ^'■'^•'(x", x'; i?), shows that the denominator for 
D > 1 does not exist. One must circumvent the arising divergence and must regularize 
properly. One chooses the coupling "7" appropriately, say, according to 

- ^ a + G(^)(a,b;0) , (2.7) 

7 

and considers the expression G^^^ (a, b; 0) — G^^^ (a, b; E) in the limit b — > a, say. We 
see that in some sense the coupling 7 has to be zero in a "suitable way" in order to 
make the final expression well defined. 

This formal reasoning can be put more rigorously [3]. As a simple example we 
start with the free Hamiltonian in IR^. It is transformed via a momentum Fourier 
transformation into a multiplication operator, so that H = h |pp/2m. A cut-off uj 
is introduced, and the coupling 7 is made explicitly dependent on the cut-off, i.e. 
7 = j{u>). All what remains is to choose 7(0;) "in such a way that the [perturbed 
operator with cut-off u] has a reasonable and nontrivial limit [H] as we remove 
the cut-off, i.e. as u> tends to infinity." [3, p. 111]. For the free particle case we chose 

2m 1 f (fip ^ m 

|p|<'^ 



—-^ = a-\ 7T- ,„ / -p-pr = a -I ttCO , (2.8) 

7(0;) (27r)3 7| |p|2 27Th' ^ ' 
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with CK e IR. We introduce the eigen-functions 0y(p) = (27r) ^/^e ^^'^ of the 
momentum operator, and the corresponding cut-off eigen-functions 0y = Xu>4'yi where 





IpI 


< w 




IpI 





xu.(p) = <!. (2.9) 

Then {E = n^\\fi^\/2m) 

1 2m 



7(0;) |p 



2 \ It^|2 IH2 



2m 1 I o 
ckH Y I ct'^p 



nM27r)3 y|p|<, -VIpP |Pp-|k| 
(a; — >• 00) 



= a + ^^V-2mE . (2.10) 

It follows that the perturbed cut-off operator H'^ (in p-space) converges in norm 
resolvent sense to a self-adjoint operator, respectively the perturbed cut-off oper- 
ator Q-^H'^T/n ^ unitary operator, hence the path integral i^^")(x", x'; T) =< 
x"| e~^^^/^ |x' > exists and is well defined, with the corresponding resolvent kernel 
given by [3, Theorem 1.1.1, p. 113] 

(x", x'; E) = G^^ (x", x'; + (x", a; E)Gf^ (a, x'; E) 



oo / / \ \ n r r-. / \ n n— 1 



X lim 

n=l 



2m 

^\^-^y' |p|2_ ikpv-y 



G'V^(x",a;E)G'r(a,x';£;) 
G'i°)(x",x';£;)+ lim s y ^ ^ ) 3 K , ^ J 



1 2m/,,, 1 



7(u;) IpP- Ikp"^^ 

= ef(x", X'; E) + »3°'(x",a;E)Gf(a.x-;E) _ ^^^^^^^ 

ra,a 

rL°L = « + ^7=2^^ , (2.11b) 

G^^\e) denotes the free particle Green function in three dimensions. Therefore 
we have obtained the three-dimensional generalization of (2.6) for V = and we 
conclude that the perturbation expansion (2.3) together with the regularization (2.10) 
is well defined and (2.11) justifies our perturbation expansion approach, respectively 
its formal series summation (2.6) generalized to the three-dimensional case. This is 
in accordance from general considerations that the perturbed Green function for a 
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(5-function perturbation located at x = a (a e ]R^,]R^) must look like (2.11) (Krein's 
formula [3], appendix A) 

Let us note that the quantity —1/A.TTa represents the scattering length of the 
perturbed operator. 

However, an analogous line of reasoning is valid for the incorporation of potentials 
F(x). Instead of performing the usual Fourier transformation into p-space, trans- 
forming the free Hamiltonian into a multiplication operator, we perform a generalized 
Fourier transformation (decomposition) in terms of its spectral expansion, taking the 
corresponding eigen- functions (pEx with energy E\ and momentum = 2mEx/h , 
yielding a generalization of (2.11), with G^'^\E) replaced by G^^\E). 

Another way to determine r^.y{E), which also includes potential problems in 
a more constructive way, exploits that the Hilbert space of the perturbed problem 
must be properly defined, i.e. one has to find the correct self-adjoint extension of the 
heuristic operator ([2, 3, 12, 19] and references therein): 

iy-y = -— A + y(x) + V(x-a)" , (2.12) 

2m 

and A denotes the Laplace-Beltrami operator in two, respectively three dimen- 
sions. Let g{r) be a solution of the corresponding minimal (reduced) radial s-wave 
Schrodinger operator (i < A < |, /3, ry G IR, < a < 2) 

^ = -£1^ + 1^^ + '+''^"° + ^W = ^^o-(IRn, (2.13) 

where the notion H denotes the incorporation of a 5- function perturbation at r = 0, 
say {H has deficiency indices (1, 1)), and V{r) E L°°(IR"'~). We then can formulate 
the following theorem: 



Theorem [3, 12] Let 



™ ^1 \ 1 

\fr in r A = 



Ff(r) = r^ Gf\r)^\ ^ ^2.14) 

A= 1 . 



All self-adjoint extensions of the operator H are given by 

H. = -f^ + f^^ + ^-+pr-' + V(r), (2.15a) 
2m dr^ 2m r 

V{H,) = 1^ e L\JR+) g,g' e AQocOR+y, ago,x = giX H^g e L2(IR+)| 

(2.15b) 

1 3 

-oo < q; < oo, -<A<-, P,r]eM, 0<a<2. (2.15c) 
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V e L°°(IR"'') real valued, and ACiod^) denotes the set of absolutely continuous 
functions on M (here = IR^). The boundary values qq^x, gi,x are defined by 



go^x = lim 



9{r) 



r— >0+ 



g{r)-go,xG^ir) 



(2.16) 



where Gf{r) denotes the asymptotic expansion of the irregular solution Gx{r) of (2.15) 
up to order r*, t < 2A — 1. 

Generally one has G^{r) = G^°^(r) + additional terms. Two special cases of G^{r) 
can be e.g. stated for the A = ^,1, i.e. for the Schrodingcr operator in two and three 
dimensions, respectively, and all other potential terms in (2.15) equal to zero, which 
will be sufficient for our purposes. Then 



G^^r) = G'r{r) 

2 

Gfir) 



2 

m 



= — T2Vr Inr , 
nn 

mri 2mr] 
1 n-r 7^r mr 



(2.17) 
(2.18) 



Take into account the (3 contributions considerably complicates the expressions 
and will not be stated here, c.f [12]. 

We consider A = {D — l)/2 and usually we identify Fx = Fd, Gx = Gd, etc. in the 
following. Note that the additional potential V{r) does not play any further role in 
the rcgiilarization. The condition ^ < A < | can be physically interpreted as follows: 
Set = A — ^, then A(A — 1) = A^^ — 1/4 and the functional weight in the radial path 
integral corresponding to the Hamiltonian (2.15) [38] is proportional to //Ar[r^]. The 
number N corresponds to the magnetic fliix inside an infinitely thin Aharonov-Bohm 
solenoid. Only fluxes less then unity are allowed and there exists exactly one bound 
state solution and the parameter of the self-adjoint extension describes the anomaly 
of the magnetic moment of the particle bound by the 5-function-like solenoid [10]. 

Summarizing, we can state that the result of the regularization procedure (2.15) 
of the formal series summation (2.6) of the path integral perturbation expansion (2.3) 
with the heuristic expression "75(x — a)" in the Lagrangian in the path integral in 
two and three dimensions yields in an analogous way as the consideration in (2.11) 
the Green function 



x(t")=x" 



n 



x(t')=x' 



l/7-lim,,j,^<,GW(x,y;f;) 



X exp 




-x^ - F(x) + 75(x - a) 



^ Jo 



exp 



h If, 



m 



dt 



- F(x) + 75(x - a) 



dt 



x(t')=x' 
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G'(^)(x",x';£;) + 



with ri^2 given by 



J- a, a 



(2.19) 



(2.20) 



and we define the heuristic point- inter action potential 7(5 (x — a) in the path integral 
by means of (2.19,2.20) with the regularization prescription (2.15). 

A convenient way to find the function g{r) is to set g{r) = Q~^{D)G\Y,q^'^'^' {r, r; E) 
= 0~^(L>)r^G[^Q(r, r; £^), where V includes all potential terms in (2.13). Here Q.{D) 
denotes the volume of the unit sphere 5"*^^"^^. 

In three dimensions it is also possible to state the regularization rule as follows [3] 



x(t")=x" 



1/ .Te'-/« / P,%xWexp 1/ 

x(t')=x' ^ 

= G'(^) (x", x'; £;) + (r 00 (i?)) "'g(^) (x", a; E)G^^) (a, x'; E) 
with T^a,l given by 



— - F(x) + 75(x - a) 



dt 



(2.21) 



rS = « + ^r-i.G(^nx",x';£;) 

or 12 



ri2 = |x"-x'|=0 



(2.22) 



x',x"=a 



provided an explicit expression for G*^^)(x",x'; E) exists. 



3. Examples 



3.1. One-dimensional examples 

For completeness we refer some one-dimensional examples. This includes the free 
particle, the harmonic oscillator and the Coulomb potential. No new material is 
presented. The results will give an easy possibility to compare the results in one, two 
and three dimensions with each other, respectively to introduce some notation. We 
also introduce the Green functions for the free particle, the harmonic oscillator, and 
the Coulomb potential. Of course, no regularization procedure is needed. 

Let us note that an implicit equation for the time-dependent propagator corre- 
sponding to (1.1) is due to Gaveau and Schulman [25]. They obtained 

t" 

K^^\x",x';T)^K^^\x",x'-T) + il [ K^^\x" ,a;t)K^^\a,x'-T - t)dt . (3.1) 

" Jt' 
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3.1.1. The free particle. Let us start with the free particle (labeled by "(0)"). The 
path integral of the D-dimensional free particle reads as (x e JR^) 



x(t")=x" 



kW(x",x';T)= J P^x(t)expl ±^dt 



x(t')=x' 



m 



D/2 



27ri hT 



(27r) 



exp 
pexp 



(x" - x')^ 



2inT 



ip- X — xj— iT 



(3.2) 



(3.3) 



The energy dependent kernel is given by 

x(t")=x' 



Gg)(x",x';£;) = i 



im 



P^x(t) exp I ^ / x^dt 



-2-1 



^ V 27r i /i 



15/2 



m 
2^' 



rr i\2 



x(t')=x' 

i(l-0/2) 



Ki D 



Ix" -x'l 



n \2n^ 
For D = 1 one has 



. / \ D/2 / / NO 

1 / m \ / m(-7r)^ „ ^ 



2E 



X — X 



U^-D/2) 



H 



(1) 

^ 2 



Ix" - x' 



V-2mE 
V-2mE 



(3.4) 
(3.5) 



From the general theory we now have [7, 27, 32, 51, 52] 

/,/f\ It 
X\t )—X 



'Dx{t) exp 



m 



+ ^5{x — a) 



dt 



x{t')=x' 



m 



hV 2E 



exp 



exp 



+ 



m7 



V—2rriE . ,, . 

z \x — X 

n 



{\x" -a\ + \a- x'\) 



2r 



(3.6) 



To determine the propagator from (3.6) the inverse Laplace-Fourier transformation 
must be applied, which can be explicitly done with result [27] 

x(t")=x' 



x'^ + ^5{x — a) 



dt 



x{t')=x' 



m 



27rinT 



exp 



im 
2hT 



/II l\2 
(X — X ) 



X erfc 



m .. . 1 

— — - a; — a\ + \x — a\ — —■yT 



(3.7) 



3.1.2. Tiie harmonic oscillator. The harmonic oscillator (labeled by "(<^)") path 
integral is given by [23] 



x{t")=x" 



x{t')=x' 



• t 

J 2n. Jfi 



mu! 



2tt i % sin a;T 



1/2 



exp 



rriijo 
2in 



(x'^ +a;"^)cota;T-2 



sin (jjT 



(3.8) 



We do not discuss the case of caustics etc. The energy dependent kernel can be 
evaluated to be given by (c.f. [5], and appendix 1) 



x[t )=X 



T>x{t) exp 

x{t')=x' 

\D_i,^ 

2 Tiuj J 2 + 



1 m 

2n Jh 



{x^ -uj^x'^)dt 



2mu 



(3.9) 



The general theory then gives for a perturbed harmonic oscillator [32] {x" > a> x'): 

x{t")=x" 



-X 



u X + 'y5(x — a) 



dt 



x(t')=x' 



m 



V2 hcv J 



2m<jO , 1 _ 
-^x' \D 1 , E 



2mijO 



h 



X 



X 



h V Trhu! V 2 hu> 



X D 



2miv 



a \D_i 



2m<jO 



a 



(3.10) 
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(x" >x'>a): 

x{t")=x" 

dT e'^^/^ / Vx{t) 



exp 



h .L, 



771 . 2 7Ti 2 2 r/ \ 

—X — —u X +^d{x — a) 



dt 



I m fl E\ f f2rruJ 



x']d 1 



2 """E^ 



2mw 



h 



X 



X 



2mu) 



a \D 



1 - ^ 



7 / m 



^ V TT^o; V 2 ^a; 



1 E 



2mu; 



h 



x"]D 1 



2mu; , 



X -D_i , _^ 

2 "I" Eo; 



2mc<; , 

a ] D 1 I B 

Tj, / 2"l"fta. 



2 mo; 



a 



(3.11) 



{x" <x' <a): 

x{t")=x" 
x{t')=x' 

m7 



TTZ . 2 2 2 r / \ 

— a; u) X -\- ^oix — a) 

2 2 ' V V 



2mc<; 



x"\D 1 



mu) 

'2n 



X 



+ 



X 



2mc^ 

— — a I — 



2muj , , 

/ 2 + hw 



2muo 



h 



X 



h V Tr/io; \2 /la; / 



1 £; \ 



- 



X D 



1,E — 



2muj 



a \D 



2muj 



n 



a 



(3.12) 



3.1.3. The Coulomb potential. The Coulomb potential (labeled by "(C)") in one 
dimension is also know as the Kratzer potential. Its corresponding Green function via 
path integration can be done by considering the £)-dimensional Coulomb potential, 
say, and then restricting to -D = 1 [17]. Of course, the problem is closely related to a 
Kustaanheimo-Stiefel transformation in the in the path integral, hence to the space- 
time transformation technique in path integrals. This has been discussed by many 
authors, let us mention [17, 21, 33, 38, 42, 44, 47, 54, 63] and references therein. A 
comprehensive survey will be given in [40]. 
We now have (x > 0, k, = qiq2-\/ —m/2E /%) 

x{t")=x" 



dTe 



iTE/n 



I 



T>x{t) exp 



x{t')-- 



m 



91 92 

\x\ 



dt 



11 



(3.13) 



For a 5-perturbed 1/r-potential we then obtain [32] (k = qxq^yj —mj^E /%, x" > a > 
x') 



X{t )=X 



n 



Vx{t) exp < — 



x{t')=x' 



^ .If 



m .2 q c-/ \ 

—X H \-jdix-a) 

2 X 



dt 



(3.14) 



and similarly as in (3.11,3.12) for the other cases. 



3.1.4. Multiple 5 -function perturbations. For completeness we also cite the case 
of a multiple 5-function perturbation. Since only the Green function G^^\E) 
of the unperturbed problem is relevant, the incorporation of multiple 5- function 
perturbations can be successively done with result [3, 28, 36] 



X[t )=X 



/ dT ^^^^1^ 
^ Jo 



Vx{t) exp ^ — 



N 



x{f)=x' 

G^^\x",x';E) G^^\x",ai;E) 
G^^\ai,x';E) G'(^)(ai, ai; E) - I/71 ... 

GW(aiv,x';E) G^^\aN,a^) 



i=i 



dt 



G^^\x",aN;E) 
G'(^)(ai,a^;E) 

G^v\aN,aN;E)-l/^N 



G'W(ai,ai;£;)-l/7i 



G^''\ai,aN;E) 



G^^\aN,ai-E) ... G^^\aN,aN-,E) - 1/^^ 



N 



= G(^\x",x';E)- Y: irf2}^:G(^\x",ar,E)G(^\ar,x'■,E) , 



(3.15) 



(3.16) 



with the matrix tJ^J^^^^ given by {{a} = {7fc}f=i, {a} = {ak}k=i) 



^J3 



(3.17) 



The determinant and matrix representations are checked by induction [36] . 



12 



3.2. Two-dimensional examples 



3.2.1. The Free Particle. We keep the ordering schema of the first subsection and 
start with the free particle case. Form the general D-dimensional free particle Green 
function we consider D = 2 and obtain 



(x",x';£;) 



Aol ^ lx - X 



n 



(3.18) 



which is logarithmically divergent for |x" — x'| ^ due to Ko{z) oc — lii{z/2) + \E'(1) 
{z 0). Here = 0.5772156649 01532 86061... denotes Euler's constant. In 

order to regularize the problem we apply the theory of section 2. We set (r > 0) 



/ N ^ rr^ (\/-2mE 
9\T) ^ -^\/rKo \ r 



from which follows 



m 



n 



^1,2 = *(!)- In 



V-2mE 



nh'' ' "^'^ '''' "' 2% 
(x e IR^ \{0}). We therefore obtain (compare also [2, 3]) 

x(t")=x" 



m 



— i +75(x-a) 



dt 



x(t')=x' 



V-2mE . „ 



+ 



m 



2^0 



V-2mE , A^^ / y/-2mE 

x — a| IAqI : la — X' 



a + 



In-^^^^ *(1) 



2h 



r 



< 



+ 



a + 



m 



In-^^ ^'(1) 

2h ^ ' 



(3.19) 



(3.20) 



(3.21) 



and polar coordinates about |x — a| have been used in the last line. The one bound 
state wave-function 



^(") = ^i^o(|x - a| e-^t^-'^V— *(!)] ) 



(3.23) 
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has the energy 



m 



(3.24) 



The normalization of the bound state wave-function foUows from the integral repre- 
sentation [29, p. 672] 



^0 



xKi,{ax)Ki,{hx)dx 



2sin(zy7r)(a2-62) ' 



> , 
zy = . 



(3.25) 



3.2.2. The harmonic oscillator. In order to discuss the harmonic oscillator we need 
the radial Green function for the harmonic oscillator in D dimensions. We have [6, 
26, 38, 55] 

f oo 



1-D 1 



^r'r")^- I dTe^^^/^ 
Jo 



r{t")=r" 

J Vr{t) exp 

r{t')=r' 



X 



1 

n 

■w 



m 



(1 + 



D-2\2 



1/4 



m 



ij?r^ 1 dt 



2mr2 



(3.26) 

We consider D = 2. For the determination of the behaviour of the Whittaker functions 
for 2; — > use has to been made of the representations [29, p. 1059] 



r(i-/.-A) 

respectively [29, p. 1063] 

(_l)2At2M+l/2g-2:/2 



(3.27) 
(3.28) 



Wx,^{z) 



r(^ + A;- A + i) 



2/ 



^ k\{2n + ky. ^ 



X 



*(/c + 1) + *(2/x + + 1) - *(/x + - A + I) - In^] 



2^ 'sr' r(2/.-/c)r(/c-/.-A + i) 



+ i-z)-^^ Yl 



(3.29) 



In the latter equation it is required that 2// e INq, and the last term is understood to 
be ignored for fj, = 0. Putting the interactions center at the coordinate origin, we set 



9{r) 



r(l - E/2hiv) 
27r?kor 



( 2X1,^ ( 

W^£/2ftu;,0 ( 1 Af£/2?iu>,0 ( 



(3.30) 
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and therefore we obtain due to the general theory 



no .2 "^22, sr \ 

— X - X +75(x) 



x(t')=x' 



r(i -^/2nw) 

2TThu>r'r" 



mu 



mo; 



+ 



WE/2huj,o\^—r jWE/2hu>,o\^j-r j 



ojr'r" 



a + 



m 



* — +ln— 2*(1) 

2 2hio h ^ ' 



1=1 



2M\T'T" ^^/^'^-.'/HT 



mu) ,2 1 , , 

r I ME/2nuj,l/2 



mu) ,,2 



(3.31) 



The bound state energy levels E^^^ are determined by 



a + 



m 
2^ 



2 2nuj 



+ ln!^-2*(l) 



(3.32) 



3.2.3. The Coulomb potential. In order to treat the two-dimensional Coulomb 
potential case, we must know the radial two-dimensional Coulomb Green's function. 
We have (c.f. the references in the one- dimensional case) 



x(t")=x" 



exp 



x(t')=x' 



^ .2 , Q1Q2 

2 X 



dt 



/ m T{1 + k) 



n 



1=0 

X 



2 W — 



(3.33) 



(k = qiq2^/—m/2E/h) and the S'|'^''(n) are the (real) hyperspherical harmonics on 
the S"*^^"^-* -sphere taken at the unit vector O on the sphere. We consider D = 2, and 
putting the interactions center at the coordinate origin, we obtain in the regularization 
procedure by means of (3.27,3.29) 



9{r) 



277% 



m r(i — K 



2E 



(3.34) 
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Therefore we get for the (5-perturbed two-dimensional Coulomb potential problem 



x(t")=x" 



^ Jo 



^T^^ET/n j p2^^^x(t)exp 



.2 . Q1Q2 ^ . 

2"^ + 



dt 



x(t')=x' 



m r(i — K 



+ 



+ 



m 



2E 2nnVr' 

2 ^2 a ^K>o(y=8S0^^«,o(V=8^^) 



2 / ^^/^// 



rr" m 
a H T, 



^{ +ln^^^^ 2*(1) 



1 



27Th V 2£; 



, -W^^ z V-SmE -^]M^i[ V-SmE ^ 



(3.35) 



The bound state energy levels E'i"^ are determined by 



a + 



m 



2T^n' 



2 n 



qiq2 / TO 



2£;i"^ 



+ ln^^ _2*(1) 



= 



(3.36) 



3.2.4. Multiple S-function perturbations. Similarly as in the one-dimensional case, we 
can consider the case of multiple (5-function perturbations in two dimensions, therefore 
a lattice of 5-function perturbations. Repeating the regularization process for each 
additional 5-function we obtain by induction [2, 3] 



x(t")=x" 



1 

^ Jo 



/ 



P^(o) x(t) exp { - 



x(t')=x' 



yx +2_^7fed(x-afe) 



N 



k=l 



m / \/—2mE 



n 



x"-x'| 1 + 



m 



V-2mE 

n 



x"-aJ \Ko 



h Jt> 

2 N _^ 
E (r{a\,{a}) . 

V-2mE . 



dt 



h 



(3.37) 



with the matrix r|f ^ {a} by [{a} = {ak}^^i, {a} = {afc}^i. Note that 

G<°'(x, y; E) = Gf\x- y,0; E) = G<°'(x - y; E)] 



m \/—2mE\ 



(3.38) 
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f + lim [G^"^ (x; 0) - G^"^ (x; E)] j = f 

1 -Gf (a,-,a;.;E) ' (3.39) 

where G^°^(x;0) = -(m/Tr/i^) In |x| (x e IR^\{0}), and &^\^;E) = G^°^(x;£;) 
for X 0, and G2°^(x;i?) = otherwise. The bound states energy levels E'n"^ are 
determined by det (rj^j^^^^^l^i"^)) = 0. 



3.3. Three-dimensional examples 

3.3.1. The free particle. Again we start with the free particle. Form (3.4) we know 
the free particle Green function and we obtain according to the general theory 



9(1-) 



m 



-r^J-lmEjh 



(3.40) 



(r = |x|, X e IR^ \{0}). Therefore we get (c.f. (2.11) and [3]) 



ir 



x(t"), 

VLo) x(t) exp 

CKjCl I ft J ^1 

x(t')=x' 

Gf (x",a;£;)Gf (a,x';E) 



-X +7()(x — a) 



dt 



a-\ — 



(3.41) 



Expanding into three-dimensional polar coordinates we obtain (compare also [3]) 



r{t") 



f/\ It 



1 ' ^ fOO 



1=0 
1 



m 



+ 75(r) 



r{t')=r' 



2iTfir'r' 



m 
2E 



sinh ^V— 2m£' exp ^ — \/—2mE 



+ 



exp ( — ^/—2mE 



r' + r' 



271;% J rr 



m I — 

a H itv —2mE 

2TTTr 



m 



h 



1,^1 ( V-2mE ^ f V-2mE ^> 



Denoting the free three-dimensional propagator by 



h 



K(0)(x",x';T) 



27rinT 



3/2 



exp 



[2hT^ ' 



(it 



(3.42) 



(3.43) 
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it is possible to apply the inverse Laplace-Fourier transformation and one obtains for 
the perturbed propagator [60] 



x(t")=x" 



Vt,o) x(t) exp <^ - 



x(t')=x' 



^ .h' 



= Kf (x",x';r) + 



X + |a-x'| + |x" -a|,0;T)dw 



— X +75(x-a) 



dt 



-2nahu/m 



{u+ |a-x'| + |x" -a|) 



'(0)^^// ^>.r^^^ K(o)(|a-x'| + |x"-a|,0;T) 



(x",x';T) + 



mla — x'l |x" — a| ^ 



= k|°)(x",x';T) + *(«)(x')*(")(x")e'^'"'^/'' 

+ 1 rrr^^ r / e'^'^l^l^'^/^fw - la - x'l - |x" - a|) 

|a-x'||x" -a| Jo 



(3.44a) 
(3.44b) 



X Kf^w- |a-x'| - |x" -a|,0;T)dw , 



(3.44c) 



for q: > 0, ct = and ct < 0, respectively. For o: < there is one bound state 
wave- function 



^(«)(x) 



Q,^2 ^-2'nah\J^-a\/m 



m |x — a| 

and the energy eigen- value i?^"^ has the form 



(3.45) 



(3.46) 



3.3.2. The harmonic oscillator. In order to discuss the harmonic oscillator we use 
the representation (A. 7). Applying the regularization procedure and by means of 
(3.27,3.29) we set 



r[|(3/2-£;/M] 



WE/2nu.,i/J-^r^]ME/2n^,i/J-^r^] ■ (3.47) 



We then obtain by denoting the three-dimensional harmonic oscillator Green function 
by G^^\E) from (A. 7) {y = — | + E/huj, a — |a|, d — y^2?mJ7^a, compare also [45]). 

x(t' 

I dT f^^T/h 

n 



' . 9 //ton ^ , X 

2"^^ - X +7<^(x-a) 



dt 



x(t')=x' 



G('")(x",x';E) + 



G'('^)(x",a;£;)G'('^)(a,x';E) 

rtJ(i) 



(3.48a) 
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where 

2(27r)W [a 



+ 



2m<jO 



n 



D'l{a)D,{-~a) + 2D'^{a)D'^{-~a) + D,{a)D'l{-~a) 



(The case a = gives:) 

m2r2[i(3/2 - £;/na;)]2-(^/'^'^-i) 



+ 



A-K^U r'r" 



X 



E/hui-^ 



2mu , 



n 



a — 



m P^V[\{Z/2-E/huj)] 



I f. r[l(/ + 3/2-£;/M] ^ ..yn(n>> 

+ ^ na;r(Z + 3/2)(rV')2 ^4^^^' (^,0) 



X ^^S/2au;,l(/+i) (^^^ J^£;/2aa;,l(«+i) 



/ muj ,,2 
r 



(3.48b) 



(3.49) 



For a = use has been made of D^(0) = ^2"" 1'^ IV{^), and D'^^z) = uD^_i{z) - 
zDi,{z)/2^ and Y[^{9^(j)) are the spherical harmonics on the sphere. The bound states 
energy levels E^^^ are determined by rj^'^j^j (E^"^ ) = 0. 



3.3.3. The Coulomb potential. Let us denote by G'^^^E) the three-dimensional 
Coulomb Green function as given by Hostler [43] {k = qiq2-\/ —m/2E /%) 



^ Jo 



x(t")=x" 



/ 



r»'^x(t) exp 



x(t')=x' 



nJt, ['2'' ^ |x| y 



mr(i — k) 

27r77.^|x" — X 

00 I 



V-SmE ^ 



n 



1=0 n=-l ' 



1 1 rmr(i + z-K) 



X 



2E (2/ + 1)! 
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(3.50) 



According to the general theory we set (c.f. (3.27,3.29) and [3]) 



9{r) = 



2E Anhr 



n 



V-8mE 



n 



(3.51) 



Therefore we obtain (compare also [3, 49]) 



x(t")=x" 



P^(c)x(t)exp<^ - 



x(t')=x' 



— -X +—— +7()(x-a) 



= G'(^)(x",x';£;) + 
where 



rS(£;) 



(3.52a) 



^rr-^/^N mT(l ~ k) , — 



27r;i"' 



X 



2W^; , (2a)M; , (2a) - W^,, i (2a)M;;\ (2a) - (2aX^ , (2a) 
(The case a = gives:) 



(3.52b) 



m r(l - k) 

\\ 

m \^r2(l-K) 



2E 47rnr'r" '^'^l /i/^'^l ^ 



+ 



1 I V-8mE - j W^^i i V-SmE — 



a + 



mqiq2 



* 1 + ^' 2 - ^' 1 - + In , ^'^^ - \ 

I 



X W^k,/+i/2(V=8^^)m,,^+i/2(^V'=8^^^ , (3.53) 

where three-dimensional polar coordinates have been used in the last line. The bound 
state energy levels £^1"-* are determined by T^^}^y{Et^) = 0. 

3.3.4. Multiple d-function perturbations. As the last example in this sequel we again 
consider the case of multiple 5-function perturbations. Similarly as before we get by 
induction [3] 



4/ dTe'^^/"^ 
^ Jo 



x(t')=x' 



P (0) x(t) exp 

{«},{a} I 11 



m 



N 



+ ^7fe(^(x- afe) 



k=l 



dt 
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N _^ 

Gi°)(x",x';£;)+ 5^ (rfj^^^^^). .G(°)(x",a,;£;)G^°)(a,,,x';£;) . (3.54) 



=1 



with the matrix tJ*^^^^^^ given by [{a} = {ak}u=i, {a} = {^k}k=i^ note 

Gf^ (x, y ; E) = Ci'^ (x y , 0; E) ^ (x - y ; E)] 



that 



(rW,{a}),,,' = + ^-^^^y^^' - Gs\^j,^j';E) (3.55) 

lim [Gf (x;0)-4°^(x;E)] j = j' , 



-Gf (a„a;.;£;) / , (3.56) 



where G^°^(x;0) = m/27rnVl (x e IR^\{0}), and Gf^(x;£;) = ^^"^(x;^) for x ^ 0, 
and G^\'x; E) = otherwise. The bound states energy levels £'n"'* are determined by 

.(0) 



«(ri"i,,.,(£r))=o. 



3.4. Nonstandard three-dimensional potentials 

In this subsection we want to discuss several potential problems important in nuclear 
physics. They are 

i) the potential well, 

ii) the Wood-Saxon potential, 

iii) and the rotating Morse oscillator. 

These potentials can serve as models for the potential strength in the vicinity of 
the nucleus more or less approximating the strong interaction force. It is assumed 
that they are radial symmetric, which allows an explicit solution for s-waves. 



3.4.1. The potential step. We first discuss the potential step. The Green function for 
the radial potential step potential V{r) = 0{b — r)Vo can be deduced from the smooth- 
step potential (see below) and has the form [37] (/c^ = 2m{E+Vo)/?i^ , = —2mE/lT?, 
I = 0) 



r(t")=r'' 



^^°°dTe'^^/^ j Vr{t)e.p ^J^^ (^^^^ -[©(,- 6) - l]Fo) c/t 



r{t')=x' 



e{b - r')e{b - r") 
Airhr'r" 



m 



2(E + Vo) 



X e 



iA:(r<-6) (^ifc(r>-fc) X + ^k 



i k{ry —b) 
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ATTfir'r" y 2E \ X-^k 

47Tnr'r" + V-E - Vq ' 

[alternatively we can write ^ — {x + ik)/x — ik) = p\. Let us denote this 
Green function by G^^^^r" ,r'; E). According to the theory we set 

a(r\- ^ I ^ p-ifc(r-b) ( Ak(r-h) ik(r-h)\(. e^ ^ ^(''-^) -p \ 

47rr/iV 2(E + ]/o) V ^ A e-^i^^'-py- 

(3.58) 

Therefore we obtain 



1 i 



r(t")=r'' 



Attt't" h Jq 



r{t')=x' 

G(^^)(r", r'; £;) + (ri7^(£;))-'G(^^)(r", 0; E)G^''^\o, r'; E) , (3.59) 



where 



+ ^-E - Vq 

and v\^^\eI^^) = determines the bound state energy levels. 



3.4.2. The Wood-Saxon potential. The Wood-Saxon potential can be seen as a smooth 
step-potential modeling the potential strength in the vicinity of a nucleus [24] . It is a 
special case of the so-called modified Posch- Teller potentials [56]. The path integral 
solution of the modified Poschl-Teller potential can be done via the path integral 
solution of the SU(1, l)-path integral [46], form which the path integral solution of the 
smooth-step potential can be derived [31]. 

From the Green function representation of the modified Poschl-Teller potential, 
respectively the Wood-Saxon potential [34, 48], 

r(t")=r" ^„ 
r{t')=r' ^ 

_ 2mR r(mi)r(mi + 1) 

~ r(mi + m2 + l)r(mi -m2 + l) 



m 9 

— r^ + - 
2 1 + e( 



-h)/R 



dt 
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X 2-Fi mi, mi + 1; mi - m2 + 1; 



1 — tanh ^2r' 



1 + tanh 

X 2-Fi mi, mi + 1; mi + m2 + 1; — — 



(3.61) 



(mi,2 = V2^R{V-E -Vo ± \/^) /n) we then can state the Green function of the 
rotating radial Wood-Saxon osciUator which has for s-waves the foUowing form [36] 



r(t )=r 

,,E.,n j 

r(t')=r' 

mR r(mi)r(mi + 1) 

~ 27Th^r'r" r(mi + ma + l)r(mi - ma + 1) 



m 



Vo 



hjf,\2 1 + e^^-b)/R 



dt 



X 



X 



1 — tanh ^^-|^ 



1 - tanh 



2 / 1 + tanh ^-|^ 



2 / 1 + tanh ^ 



2i? 



2 y V 2 

X { 2F1 ( mi, mi + l;mi - m2 + 1; 





2 




+m2 




2 




-6 X 



2R 



l + tauh 

X 2F1 mi, mi + 1; mi + m2 + 1; r — — 



1 — tanh ^ 

2-ri mi, mi + 1: m.i + m.2 + 1; — 



1 + tanh ^ 
2-Fi mi, mi + 1; mi - m2 + 1; — — 



1 — tanh 

X 2F1 { mi, mi + l;mi - m2 + 1; — — 



1 — tanh ~° 

X 2F1 ( mi,mi + l;mi - m2 + 1; — ^ 



(3.62) 



Let us denote this Green function by G'^^^^ {r" ,r' \E'). Following the general theory 
we set g{r) = rGl^^\r, r; E). Let us abbreviate 



, litaiili^TT^ 
f±Kf) = 2F1 mi, mi + 1; mi ± m2 + 1; — 



(3.63) 
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Note the relations 



(1 T tanh — \ 
mi,mi + l;mi±m2 + l; ^^^) ' 

f _ j_2Fi{mi,mi + l;mi ±m2 + 1; 1(1 ^tanh 3^)) 

4i?cosh2^ ' 
_ 2F{'{mi,mi + I;mi±m2 + 1; |(l^taiih^)) 
•^^^ ^ ~ 16i?2cosh^^ 



(3.64) 
(3.65) 



16i?2 cosh* 2^ 
6 2-P'i'("7,i,mi + l;mi ±m2 + 1; Kl^Ftanh^^)) 

tj&illll — ~r o — ; 

2i? 4i?2cosh^2^ 



(3.66) 



We now find 



yo,3 



i?r(mi)r(mi + 1) 



r(mi + m2 + l)r(mi — 7712 + 1) 



X 



1 + tanh 2^ \ "^^-'"^ / 1 - tanh ^ \ '"^+™^ 



/_(0)/;(0)-/i(0)/+(0)J 
(3.67) 



Aws) 

91,3 



RT{mi)T{mi + 1) 



r(mi + m2 + l)r(mi - m2 + 1) 



X 



X 



1 + tanh ^ \ "^^-^^ / 1 - tanh ^ \ "^^+^^ 



^(/-(o)/^(o) - r (o)/+(o)) + /:(o)/;(o) - ^/-'(o) 



(3.68) 



Therefore we obtain for the 5-perturbed Wood-Saxon potential 



4:7rhr 



r(t J=r 

l;:i?jfdTe^^^/^ [ P,c..)r(t)exp 



'•{t')=r' 



m 



Vo 



2 + 1 + e(^-b)/^ 



dt 



= G^'^^\r",r';E)+{T%'\E))-'G^'^^\r",0;E)G^'^^\0,r';E) , 



(3.69) 



with T%^\e) given by 



i a,0 l-^) - «S'0,3 - 9l,3 , 

and the bound-state energy levels E^ are determined by T^^Q^\Et^) = 0. 



(VKS) {WS) 



(3.70) 
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3.4.3. The rotating Morse oscillator. Another model for a nuclear potential is the 
so-called rotating Morse oscillator [59]. Its Green function can be derived by path 
integration via means of the Green function of the Morse potential [14, 20, 30, 54] and 
has for s-waves the following form [36] 



ATrhr'r" 



Vr{t) exp { ^ 



m .2 FoV 

— r 

2 2m 



(e-^'"-2ae-'^) 
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m 



rU + V-2mE/h-aVo) 



AnVoh^r'r" r(l + ^/-8mE/h) e(^'+'^")/2 

X I ^aV^o,V^/R(2^oe-^< )M^^^_^/^/,(2Vbe-'"> ) 



W. 



aVo,V-'2rnE/h 
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^aVo,V-'2mE/h 



(3.71) 



Let us denote this Green function by GQ^\r" , r'; E) 



Mr) ^ M^y^^^^/^{2Voe-n , Mr) ^ W^y^^^^/^{2Voe-n • (3.72) 



Note the relations 



MO) = M,y^^^^/f,{2Vo) , ^(0) = -2VoM;^^^^3^/,(2Vo) 
w'liO) = 2VoM'^y^^^^^^{2Vo)+W,'M'J^y^^^^^^{2Vo) . 



(3.73) 
(3.74) 



According to the theory we set g{r) = rG^^^r, r; E). Therefore we obtain 
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l(^i(0)^'(0) -<(0)«;2(0)) +^(0)^(0) - ^^'(0) 



(3.76) 
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Therefore we obtain for the (5-perturbed rotating Morse oscillator 



AttUt't" Jq 



/oo 




2a e"^) dt 



r{t')=r' 



= G'(^) (r", r'; E) + {T%^ {E))~'g^^^ {r", 0; E)G^^^ (0, r'; E) , 



(3.77) 




(3.78) 



and the bound-state energy levels En are determined by F^^^q {En ) — 0. 
4. Summary and discussion 

In this paper I have presented a path integral approach for the incorporation of two- 
and three dimensional 5-function perturbations by means of a perturbation expansion. 
The motivation was to demonstrate the validity and applicability of the formalism in 
the context of path integrals, i.e. "to build up quantum mechanics from the point 
of view of fluctuating paths [21]", hence to put the results of the operator and the 
path integral approach in quantum mechanics on an equal footing, and to supply 
steps towards a complete classification of solvable Feynman path integrals [39]. A 
regularization procedure was needed to treat appearing divergencies, but could be 
systematically discussed. As was shown, a regularization procedure known from 
functional analysis could be applied in the context of path integrals. In order to do 
this, the coupling 7 in the heuristic expression "75(x — a)" in the Lagrangian in the 
path integral had to be properly interpreted and it was shown that the formal series 
summation (2.6) of the perturbative path integral approach (2.3) corresponding to 
the Hamiltonian (2.15) could be regularized for two and three dimensional 5-function 
perturbations, making the two and three dimensional analogue of (2.6) well-defined. 
We succeeded in developing the theory intaking into account potential problems 
together with 5-f unction perturbations in two and three dimensions. The general 
feature of the Green function of the perturbed problem with one point-interaction at 
X = a (x', x", a e H^, H^) has the form 



G(^)(x^x';£;) = G(^)(x^x';£;) + (rOO(£;))-'G(^)(x'',a;£;)G(^)(a,x';£;) (4.1) 



giving in general a transcendental expression for the En . Generally, a <0 generates 
bound states, whereas a > generates resonance states. 
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We demonstrated the formalism with a sample of two- and three-dimensional 
examples. They included the free particle, the harmonic oscillator, the Coulomb 
potential case, and some potential models familiar in nuclear and electronic shell 
physics. The latter were studied in three-dimensional space for s-waves only. We 
chose the potential step, the Wood-Saxon potential, and the rotating Morse oscillator. 
Of course, the Wood-Saxon potential is just but one model in the class of potentials 
related to the modified Poschl- Teller potentials, other exactly solvable potentials can 
also be considered, e.g. the modified Poschl- Teller potential itself, or more specifically 
the Hulthen or a Scarf-like potential. These latter potentials can be used to describe 
a screened Coulomb potential, where an additional 5-function perturbation at the 
center would model a point interaction (of the nucleus) in the environment of an 
overall Coulomb-screening, say, of an electron shell, respectively, a point interaction 
(inside a nucleus) in the environment of an overall screening of the strong force in the 
nucleus itself, and many more. 

One may ask, what happens in higher dimensions? Actually, point interactions 
do not make sense for D > 4: "Any possible mathematical definition of a self-adjoint 
operator H of the heuristic form —A + X5y in L^(IR'^) should take into account the fact 
that, on the space Cq°(JEI'^ \{0}) of smooth functions which vanish outside a compact 
subset on the complement of {y} in IR^, H should coincide with —A. For d > 4 this 
already forces H to be equal to —A on if^'^(lR'^) since — A|(-;oo(-jf^d y^g}) essentially 
self-adjoint for d> 4 [3, p. 2]", and c.f. [58, Theorem X.ll, p. 161]. What remains are 
5- function perturbations on planes and hyperplanes, respectively, along perpendicular 
lines and planes [36], systems which can model a Casimir effect, e.g. [11]. 

In one dimension, making the strength of the 5-function perturbation infinitely 
repulsive, produces Dirichlet boundary-conditions at the location of the 5-function 
perturbation. In two and three dimensions, in setting a — some sort of defect 
is produced, which one may call a zero-radius hard core. Let us emphasize that 
again the entire Green function has to be taken into account and not only part of it 
(compare [50]). This shows once more that it is not enough just to throw away the 
continuous part of the spectrum (say, in the Coulomb case), but a proper regularization 
always takes into account the entire Green function yielding inadmissible contributions 
from both the discrete and continuous spectrum of the unperturbed problem. This 
concludes the discussion. 
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Appendix 1. Green function for the harmonic oscillator 



We want to demonstrate how to calculate the Green function of the harmonic 
oscillator. Of course, it is always possible to determine the Green function of a 
quantum mechanical problem by means of an operator approach exploiting the theory 
of differential equations. However, this is not the point of view here. The purpose here 
is to show how a Green function can be obtained by knowing the propagator form a 
path integral calculation. 

We consider the propagator K{x" , x'; T) (3.8) of the harmonic oscillator. We need 
the integral representation ([13, p. 86; 29, p. 729], ai > a2, 3?(| + /x — i^) > 0): 
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Here Wi,^fj,{z) and M^^^{z) denote Whittaker-functions. Furthermore we make use of 
the relations of the parabolic cylinder functions in terms of Whittaker-functions [13, 
pp.39] 
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We now obtain 
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(perform coordinate transformation and Wick-rotation) 



I J.I I roo / 

-Jo 



2%' 



coth 



2 J 



28 



X -'i/2l-^a;a; smnv 1 + /_i/2 I -^a^ smnv 
(apply (Al.l), set v = -\^ E/Tiw) 

1 f 2 



X 



(apply (A1.2)) 
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X 



(A1.3) 



which is the result of (3.9). The particular structure of the propagator allows to 
derive the higher dimensional case from the lower one. Let x G IR"^ and define 

2 2 

11 = -x! + x" , v = • x" , then the following relation for the Feynman kernel for the 
harmonic oscillator is valid 



K(^)(x",x';T) = ^Ak(^-2)(x",x';T) 
271 ov 



(A1.4) 



(A1.5) 



Introducing ^ = i(|x' + x"| + |x" — x'|), ry = |(|x' + x"| — |x" — x'|) we obtain for 
dimensions -D = 1, 3, 5 . . . for the Green function 
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Note the various sign conventions used in the literature [5]. 
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Appendix 2. Propagator in two dimensions 

We consider (3.21). The first term is just the free particle propagator in two 
dimensions. To treat the second term we consider the following inverse Laplace 
transformation table [22] : 
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Obviously, we will get a rather complicated expression, so that we keep the discussion 
short. We use imaginary time r = iT, which has the consequence that we must set 
—E^s. Furthermore we introduce the abbreviations 



a = 



'2m 

IT 



x" - al 



= — z — X — a 
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(3=. ^j^^-c^f^Vm-ni) (A2.1) 



Applying successively the above rules we obtain for the propagator of the free particle 
perturbed by a point interaction in (compare also [1]) 
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